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Abstract 

We discuss the non-thermal phase structure of a chirally invariant Higgs-Yukawa 
model on the lattice in the presence of a higher dimensional <h®-term. For the 
exploration of the phase diagram we use analytical, lattice perturbative calculations 
of the constraint effectice potential as well as numerical simulations. We also present 
hrst results of the effects of the $®-term on the lower Higgs boson mass bounds. 
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1. Introduction 


In this letter we investigate the influence of the addition of a diniension-6 operator 
to a chirally invariant Higgs-Yukawa model. This model can be understood as a 
limit of the standard model (SM) without gauge fields. In particular, we consider 
a complex scalar doublet and one doublet of mass-degenerate quarks. Our aim is 
to explore, whether a dimension-6 operator, for which we will employ a 
term with a coupling constant Ae, can modify the phase structure of the Higgs- 
Yukawa sector of the SM and may alter the lower Higgs boson mass bound as 
already observed in 0 , 3 . For a phenomenological analysis of a (y9'l(^)^-term see 
e.g. 0 , 3 . 

The motivation for adding a ((^l(^)^-term is twofold. First, since the Higgs- 
Yukawa sector of the SM is trivial, the cut-off cannot be removed and hence such 
a term is in principle allowed. In addition, if small values of the cut-off of 0(1) — 
—O(10)TeV are considered as done in this work, such a term can have a significant 
effect. Second, the apperance of a ((/?l 9 ?)^-term can be understood to arise from an 
extension of the SM. Studying the system with such a term could hence provide 
bounds on the couplings of such extensions in case the lower Higgs boson mass 
bound is incompatible with the Higgs boson mass of about 126GeV. The effects of 
higher dimensional operators on the vacuum stability is discussed in 0, S 3 . 

We use a lattice regularization of the Higgs-Yukawa model which eventually also 
allows non-perturbative numerical simulations for large values of Ag. The notion of 
^ exact lattice chiral symmetry 0 which derives from the Ginsparg-Wilson relation 
allows us to emulate the continuum Higgs-Yukawa sector of the standard model 


on a discrete Euclidean space-time lattice. To this end, the overlap operator [10|, lU 
as a local l3| lattice Dirac operator has been employed to study the phase structure 


of the lattice theory [13|, [1^, to derive lower and upper Higgs boson mass bounds 
M 0 0, [l0 and to analyze the Higgs boson resonance non-perturbatively (l^ . 
For a review, see j^. 

For our investigations we perform analytical calculations of the phase structure 
of the model by computing the constraint effective potential (GEP) to the first 
non-trivial order in lattice perturbation theory. In this calculation, we employ the 
same chirally invariant lattice formulation of the Higgs-Yukawa model as it is used for 
the numerical computations. We compare results for the phase structure obtained 
from numerical simulations to our perturbative predictions. In addition, we will 
provide first results for the lower Higgs boson mass bounds in the presence of the 
dimension-6 operator as obtained from the analytical, perturbative calculations of 
the GEP. 


2. Basic definitions 

In this work, we restrict ourselves to the case of one fermion doublet "0 = {t, b)'^ 
with mass degenerate quarks. The scalar fields are a complex doublet (p. Here, we 


will only provide the basic definitions of the model and refer to ref. [2^ for a more 


2 

















detailed explanantion. In Euclidean space time the continuum action is given by: 


+ J {i$t+ b(^b + y + h.c.} , (1) 

with = iT 2 (p* and T 2 being the second Pauli matrix. Besides the standard bare 
parameters and A for the Higgs potential and y for the Yukawa coupling, we add 
the dimension-6 operator Ag to the action. 

For the numerical implementation of this model we use a polynomial hybrid 
Monte Carlo algorithm (^| with dynamical overlap fermions, see ref. for details. 
On the lattice, it is convenient to rewrite the bosonic part of the action in the 
following wajlil: 

sbI*) = -k y + '*’-<‘1+E (“P'+^ -1]'+[K'-tj’). 

X 

( 2 ) 

Here the scalar held, $, is represented as a real four-vector and the relation to the 
continuum notation is given by: 


ip - \/ 2 k ( ^0 _ ^^3 


1 - 2A - Sk 


mo = 


A = 
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A 

4^’ 


Afi — 


■^6 


( 3 ) 


As said above, our main goal is the exploration of the phase structure of the model 
in the presence of the term with coupling strength Ag. We will use the mag¬ 

netization m as the order parameteJl. The magnetization is given by the modulus 
of the average scalar held and is related to the vacuum expectation value {vev) via: 


m = 




E*- 


vev 


= ■ 


m. 


( 4 ) 


For a determination and detailed discussion of the phase structure of the model 
for Ag = 0, we refer to refs. [l3|, [14 . 


3. The constraint effective potential 

Before resorting to numerical simulations, we study the phase structure analyt¬ 
ically in lattice perturbation theory for which we employ the CEP H0- We 
assume the scalar held to be in the broken phase, so the scalar held decomposes 
into the Higgs mode, h, and the three Goldstone modes, with 0 = 1,2, 3. The 


^The lattice spacing is set to one throughout this paper. 

^Here we are only interested in transitions between the symmetric and the spontaneously broken 
phases and thus will not consider the staggered magnetization 


13|,ll4|. 
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CEP U{v) is described by the zero mode of the Higgs field, Iiq = The per¬ 

turbative calculations are done by keeping the lattice regularization explicitly, i.e. 
the overlap operator is used for the fermionic contribution and all sums over lattice 
momenta are performed numerically. 

To obtain the potential the bosonic non-zero modes are integrated out. To do so, 
the bosonic action is separated into a Gaussian contribution which can be integrated 
out leading to the bosonic propagators. The remaining terms are treated as an 
interaction part and can be expanded in powers of the couplings. This separation 
into a Gaussian and an interaction part however is not unique and we employ two 
versions of the GEP. 

A derivation of such a lattice constrained effective potential can be found in 


l4 |23| . Following the procedure in these references, the Gaussian contribution to 
the action reads: 


sr'^"[h, £/"] = 2 ( ~^-P K + YI 3-p {f + Tp 

p^O V a 

which leads to the propagator sums: 


V ^ ip- -|- nP 

p^O ^ 0 


( 5 ) 


( 6 ) 


As in 1 ^, tuq is replaced by the renormalized masses in the propagator sums. The 
the mass of the Goldstone boson is set explicitly to zero. This leads to: 


1 




1 


2 ’ 


V ^ + niTT 




( 7 ) 


p^O 


The determinant from integrating out eq. 
be neglected for the GEP. 

The GEP up to the hrst order in A and Ae is then given by: 


is independent of h and can therefore 


Ui{v) = Uf{v) + Ah"^ -f Agh® 

+ A • ■ 6{Ph + Pg) + Ag • {v^ ■ (45P2 + 5APgPh + dSP^) + . (i5P^ + gp^)) . 

( 8 ) 

The fermionic contribution, Uf, originates from integrating out the fermions in the 
background of a constant held. It takes the form, 

2 


Ufiv) = log 


V 


z/+(p) + y-v- 1 - 


2p 


(9) 


where Z 2 =*=(p) denotes the eigenvalues of the overlap operator. 


= p 1 + 


±i + rfp — p 

\Jp^ + - pY 


p^ = 




sm 
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p^ = Y ■ 

( 10 ) 
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In this equation r denotes the Wilson parameter and p (0 < p < 2r) is a free param¬ 
eter of the overlap operator which can be tuned to optimize its locality properties 
12| . Throughout this work, we set r = 1 and p = 1- 

In addition to the procedure leading to Ui{v), eq. (IH]), another ansatz in per¬ 
forming the Gaussian integral is to collect all the terms that are quadratic in the 
bosonic non-zero modes from the self interaction: 


^gaussj-^^ ^ ^ ^p2 ^2 12SOAeU'^) hp 


+ 5^P-p(p' + m^ + 4At)^ + 6A6t)^)p“ . (11) 


In this approach the bosonic determinant can no longer be neglected for in potential 
calculation, since it depends explicitly on the zero mode. Further, at first order 
in A and Ae of perturbation theory, the propagator sums and combinatorial factors 
change. 


m, 


U^iv) = Uf{v) + + \v^ + 


21 / 


.2 


4x3 


log (p^ -I- mp -I- 12Af)^ -I- SOAg'O'^) • (p^ + ml + 4At)^ -|- bAg'O'^) 


-l- A I 3 Pfj 6 PfjPc T 15 Pq j -|- Agh I 45 P^j P 54 PjjPq P 45 Pq 


Ag [IhPfjP 27 PfjPc + 45 PhPg + 105 P, 


G ) 


( 12 ) 


with the propagator sums given by: 


Ph = 




V ^ p^ -I- mn -I- 12v‘^\ P SOh^Ag ’ 

p/o 


Pg = 




V ^ P mlP 4t)2A -1- Qv^Xq ' 

p/o 


(13) 

In this approach logarithmic terms appear. Depending on the choice of the 
bare parameters (mg. A, Ag), the arguments of the logarithms may become negative, 
leading to the well known problem that the effective potential becomes complex . 
We remind, that the lattice spacing is set to one implicitly such that, even though 
we use the continuum notation, all quantities are dimensionless. 

Using the analytical form of the CEP, the vev can be obtained by the (absolute) 
minimum of the potential. In order to introduce a physical scale, we set the lattice 
vev to the phenomenologically known value of 246 GeV and define the cutoff, A as 
the inverse lattice spacing: 


dP(f)) 


dh 


= 0 , 


A = 


246 GeV 


vev 


(14) 
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Further, the squared Higgs boson mass m‘jj is determined by the second derivative 
of the potential at its minimum, 


d^Ujv) 

dh^ 



v=vev 


(16) 


Due to the explicit appearance of the Higgs boson mass in the propagator sum 
eq. (I7j) for the potential Ui, eq. (IH]), we have to use an iterative approach in the 
determination of a solution for the minimum of the CEP and the Higgs boson mass. 
To this end, we hx the parameters m^, ?/, A and Ag, guess an initial Higgs boson 
mass and simply iterate eqs. 01411151) until we hnd convergence. 

We will compare results obtained from both forms of the potential to results from 
our non-perturbative simulations. As we will see below, we indeed hnd parameter 
sets, where the perturbative CEP describes the non-perturbative data well, even 
on a quantitative level. This will allow us to obtain results for the phase structure 
of the Higgs-Yukawa model considered here from the analytical perturbative CEP, 
where a non-perturbative simulation is not feasible anymore, i.e. for large lattices 
or large cut-offs. 


4. Results 

For our study of the phase structure we performed simulations for two values of 
Ag (0.001 and 0.1). Note, that having set the lattice spacing to one, Ag is treated as 
a dimensionless coupling constant. For each value of Ag we choose a set of values for 
the quartic coupling, A. The Yukawa coupling, y, is chosen such that the quarks in 
our model have a mass of that of the physical top quark, mt = y- vev ■ A 175 GeV. 
The phase transition between the symmetric and spontaneously broken phases is 
probed by scanning in the hopping parameter, k. 

In hgure [T] we show results for the bare vev computed on lattices with volume 
16^ X 32 for Ag = 0.001 (left) and Ag = 0.1 (right). Onr data show the same 
qualitative behaviour for both values of Ag. The phase transition is of second order 
when A is chosen negative and its absolute value is small. Increasing the absolute 
value of A will hnally result in a change to a hrst order phase transition. The 
appearance of these hrst order phase transitions is a natnral consequence of adding 
the dimenson-6 opeartor, (</? V) ) which can lead to multiple minima of the potential 
with non-vanishing vev. 

For Ag = 0.001 which is shown in hg. [Tal the simulation data and the analytical 
results from both versions of the effective potential agree quite well. The results 
from U 2 , eq. (El), actually coincide with the simulation data on a quantitative level 
as long as the transition is of second order. The effective potential Ui reproduces the 
behaviour of the simulation data qualitatively. However, the exact numerical results 
for the vev differ and the phase transitions are shifted to larger absolute values of 
A. 

For Ag = 0.1 which is shown in hg. Ilbl the effective potential Ui shows qualitative 
agreement with the simulations. The effective potential U 2 fails to describe the 
numerical data and the 1-loop evaluation of the CEP seems not to be sufficient. 
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(a) As = 0.001 


(b) As = 0.1 


Figure 1: Data which were obtained from numerical simulations and the perturbative approaches 
described in section |3] are compared. The plots show the vev as a function of k while Ag is kept 
fixed to Ag = 0.001 (left) and Ag = 0.1 (right) for various A. The simulation data are depicted by 
the open squares, the crosses indicate the vev obtained from C/i, eq. (HI), while the dots show the 
corresponding results from 1 / 2 , eq. (fT^ . All data have been obtained on 16^ x 32 lattices. 


The results discussed above are obtained on a relatively small lattice of size 
16^ X 32. To verify the order of the phase transitions, simulations and analytical 
calculations on significantly larger lattices are necessary. In hgure [2], we show results 
for the vev as a function of k on various volumes. The parameters are chosen in a 
region where the small volume data indicate a second order transition, fig. [211 and 
a first order transition, fig. [20 In addition, we compare the simulation data to the 
analytical results from eq. ([E]) and eq. (Ha. 

As it is shown in fig. |2al the larger volume data confirm the second order nature 
of the phase transition. Furthermore, the finite volume dependence of the second 
order transition is very well described by both versions of the effective potential. 

In fig. we show the vev obtained from the effective potential U 2 and from 
our non-perturbative lattice simulations on various volumes. Both methods give 
compatible results on a qualitative level and just the exact position of the phase 
transition is slightly altered. The jump in the vev indicates strongly the existence 
of a first order phase transition at a Ktrans ~ 0.12277. For these parameter choices, 
hnite size effects are very small. In particular, for k < Ktrans the vev stays non¬ 
zero. This means that the first order transition occurs between two minima of the 
potential with non-zero vev. Hence, this transition must occur between two broken 
phases. 

Close to the point where k ~ ^trans, tunneling events occur between the two 
minima in the simulations and hence the lattice simulation data may not agree with 
the results from the effective potential. This stems from the fact that the CEP gives 
only solutions at one of the minima and thus cannot take into account tunneling 
effects. In fig. [3a] we show the Monte Carlo time history of v at different values of 
K which clearly shows tunneling events. While for k, = 0.11757 and k = 0.11763 v 
fluctuates around the mean value of vev ~ 0.15 and vev ~ 0.40, respectively, for 
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(a) A = -0.0085 (b) A = -0.0088 

Figure 2: We show the finite volume effects of the phase structure scan for Ae = 0.001. The plot 
on the left hand side shows data for A = —0.0085 where the simulations (open boxes) indicate a 
second order phase transition. The plot on the right hand side shows results for A = —0.0088, 
where the transition is first order. In addition to the simulation data we show the data obtained 
from U 2 eq. (IT^ (dots) for both and from Ux eq. dH) (crosses) for the left plot. 


K = 0.11760 tunneling events between these two values appear, typical for a first 
order phase transition. 

From the histogram of v with an appropriate binning size, we can construct 
an effective potential from the simulation data. This is shown in £g. |3bl It is 
demonstrated nicely how the absolute minimum at around h 0.15 abruptly jumps 
to h ~ 0.35. Such a behaviour is typical for a first order transition. 





^conf 

(a) trajectories 


V 

(b) CEP from simulations 


Figure 3: The left plot shows the trajectories for ensembles generated around the first order phase 
transition generated on 16 x 32 lattices. The data correspond to Ag = 0.1 and A = —0.38. The right 
plot shows the corresponding CEP as it was obtained by taking the logarithm of the histograms 
of the magnetization. The lines in (b) just serve to guide the eye. 


Given the fact that for small values of Ae the effective potentials describe the 
simulation data on a quantitative level, it can be utilized to investigate the behaviour 
of the vev further. Due to the wider range of applicability we restrict ourselves in 














Figure 4: Here we show the CEP Ui, eq. m, for fixed Xq = 0.001 and various k values around 
the phase transition. The left plot (A = —0.0088) shows a second order phase transition for 
K 0.122715. Note that the effective potential at /t « 0.122764 actually corresponds to a crossover 
transition, see the discussion in the text and fig. The right hand plot (A = —0.0089) also shows 
a second order transition at k ~ 0.12271 and a first order transition k ~ 0.1227565 


the following discussion to the potential Ui, eq. (|H])- 

We plot the behaviour of the effective potential as a function of k in hgure|4]for a 
hxed value of Ag = 0.001. In hg. HaJ the behaviour of the effective potential shows a 
second order phase transition: the minimum moves from a zero to a non-zero value 
in a smooth way, indicating the second order nature of the transition. 

However, when A is slightly changed to A = —0.0089 we observe, in addition to a 
second order transition at k ~ 0.12271, a phase transition from one non-zero value 
of the vev to another non-zero value of the vev at large k- values, as shown in £g. |4bl 
This transition happens through a double well potential which is almost realized at 
K = 0.1227565. 

To determine the location of a second order transition in the CEP, we investigate 
the curvature of the potential at its minimum, U {vev). The curvature of the 
potential in its minimum is related to the susceptibility x of the magnetization, y oc 
1/U {vev), and is therefore minimal at the location of the second order transition. 
The susceptibility at the phase transition diverges when the volume goes to inhnity 
corresponding to U {vev) going to zero. To study this hnite size effect, we investigate 
the behaviour of the vev and the inverse curvature of the potential for volumes up 
to 128^ X 256. Some example plots are shown in £g. |5] where we plot 1/U {vev) as a 
measure of the magnetic susceptibility. In hg. EUthe typical behaviour for a second 
order transition is apparent for A = —0.007. For A = —0.0085 (hg. l5b|l a second 
maximum in the inverse curvature of the potential is visible. This second maximum 
is volume independent and indicates a crossover transition in the broken phase. In 
hg. IS^the second transition aX k = 0.12275 has turned into a hrst order one, while 
the second order transition between the symmetric and broken phase is still present 
at smaller values of k. 

Our results for the phase structure computed within the framework of the CEP 
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Hi 


(a) A = -0.007 (b) A = -0.0085 (c) A = -0.009 


Figure 5: Here the volume dependence of the location of the minimum of the CEP Ui, i.e. the 
vev (upper plots) and its inverse curvature in the minimum as a measurement for the magnetic 
susceptibility (lower plots) are shown as a function of k for Ag = 0.001 and a set of A-values. 


are summarized in fig. IH] for both Ag values. For Ag = 0.001 we clearly observe a 
second order phase transition at small absolute values of A. At intermediate absolute 
values of A an additional crossover transition sets in within the broken phase. This 
crossover turns into a hrst order phase transition around A ~ —0.0089. The second 
order transition still exists at this point separating the broken and symmetric phases. 
Around A ~ —0.0098 and k ^ 0.12267 the line of second order transition runs into 
the line of first order transition. From that point on only the first order transition 
remains separating the symmetric and broken phases. 




(a) Ag = 0.001 (b) Ag = 0.1 

Figure 6: Phase structure obtained from the CEP C/i (l8|). There are two phases - a broken and a 
symmetric one - separated by lines of first and second order phase transitions. Furthermore there 
is a small region in parameter space, where a first order transition between two broken phases 
exists for Ag = 0.001 and Ag = 0.1. The lines between the data points are just to guide the eye. 

For Ag = 0.1 the general behaviour is very similar. However, the region in 
parameter space where the additional transitions between two broken phases occur 
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is extremely narrow, see the inlet in fig. [60 In fact, the region is so narrow that it is 
well possible that in infinite volume only a single transition line exists with second 
order transitions for larger and first order transitions for smaller quartic couplings. 

With the CEP the Higgs boson mass can also be obtained from eq. (lT5l) . In 
figure El we show some first results for the cut-off dependence of the Higgs boson 
mass obtained by the CEP Ui for a series of A values around the region where 
the first order transitions appear. For Ae = 0.001 we observe, see fig. (70 that for 
the range of cut-off values considered here, the Higgs boson mass can be lowered 
compared to the lower Higgs boson mass for vanishing self couplings A and Ae as 
was also found in ref. [ij. 

Inspecting, however, fig. [70 we find that for Ag = 0.1 and for small cut-off values, 
the Higgs boson mass is significantly larger than the lower bound at vanishing A and 
Ag. Note that mn/A 0.1, i.e. we are still staying in the scaling region of the model. 
The increase of the Higgs boson mass at small cut-off can be understood from the fact 
that the Ag(<h^<h)^ term in the action provides a positive contribution to the Higgs 
boson mass shift, dominating the negative contribution from the Yukawa coupling. 
For larger values of the cut-off, the Ag coupling becomes less and less relevant and 
the Yukawa term provides the major contribution to the mass-shift such that we 
eventually find the standard behaviour of the Higgs boson mass as a function of the 
cut-off in fig. [70 

We plan to investigate the cut-off dependence of the Higgs boson mass through 
non-perturbative numerical simulations in the future. However, if the picture of 
fig. [70 is confirmed, this would lead to a bound on the values of Ag since the 126 GeV 
Higgs boson mass would be in conflict with the cut-off dependent mass at low values 
of the cut-off. As a consequence, only rather small values of Ag oc 0(0.001) would 
be compatible with the 126 GeV Higgs boson mass. 
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Figure 7: Shown is the cut-off dependence of the Higgs boson mass obtained from the CEP 
according to eq. (1T71) for A = 0.001 on a 64^ x 128-lattice (left) and A = 0.1 on a 192^ x 384 (right). 
In both plots we also show the standard model lower mass bound (Ag = A = 0). 
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5. Conclusions 

In this letter we focused on the investigation of the phase structure of a chirally 
invariant lattice Higgs-Yukawa model including an additional higher dimensional 
operator, with coupling strength Ag in the action. For the analysis of such 

a system we restricted ourselves to small values of Ag for now. This allowed us 
to compare our numerically obtained results with analytical predictions from the 
constraint effective potential evaluated in the same lattice setup as the numerical 
simulations were carried through. 

In general, we obtained a very good qualitative and even quantitative agreement 
between both approaches leading to the phase structure shwon in £g. [U] for hxed 
values of Ag = 0.001 and Ag = 0.1. 

Fixing Ag > 0 stabilizes the potential, allowing thus to drive the values of A more 
and more negative. For sufficiently small values of A we observe smooth transitions in 
the magnetization, fully compatible with the second order phase transitions observed 
for Ag = 0. However, from a certain negative value of A on, we hnd an additional 
phase transition which can be a crossover or first order transition. Indications for 
these transitions can be detected from the behaviour of the magnetization computed 
both in the effective potential and the numerical simulations, see e.g. hg. |4bl Thus, 
the resulting phase diagram in fig. [6] turned out to be rather rich with second and 
hrst order phase transition lines when changing k. We note in passing that by hxing 
the hopping parameter k and hence the bare Higgs boson mass, it is possible, to 
move to a broken phase by only changing the quartic coupling of the theory. 

A natural extension of the investigation here would be the exploration of the 
phase structure of the model at non-zero temperature. Our results show that a 
simple extension of the Higgs-Yukawa sector of the standard model by a (0V)^ 
term leads to hrst order phase transitions. This might open the possibility to gen¬ 
erate a strong enough hrst order pha se transition at a non-zero temperature which 
is compatible with baryogenesis [2^ even at a value of the Higgs boson mass of 
126 GeV. 

The constraint ehective potential also allows to compute the Higgs boson mass 
from the second derivative at its minimum. By hxing the value of Ag = 0.001 and 
driving A more and more negative, we obtain lower and lower values of the Higgs 
boson mass and, in particular, substantially smaller values than obtained for Ag = 0 
at a comparable value of the cut-oh. This hnding is fully compatible with the results 
of [l|. As a criterion to obtain an absolute lower bound for the Higgs boson mass 
one may choose the value of the quartic coupling, where the second order standard 
model like phase transition turns into a hrst order one since in the Higgs-Yukawa 
sector of the SM itself only second order phase transitions occur. 

We have also found that for larger values of Ag = 0.1 and at small values of the 
cut-oh the positive contribution of the Ag term to the Higgs boson mass-shift leads to 
signihcantly enhanced Higgs boson masses. In fact, we can already exclude certain 
values of the quartic and Ag couplings since there the 126GeV Higgs boson mass is 
in conhict with the lower bounds obtained here. It will be interesting to perform a 
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more systematic study of the lower Higgs boson mass bounds at additional values of 
Ag. By employing also numerical simulations this can provide exclusion bounds for 
the coupling values and hence for models which lead to an extension of the standard 
model with a term. We plan to carry out such investigations in the future. 
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